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^ ■ Abstract 

£f) ■ We have carried out a systematic analysis of the transverse dipole spin re- 

sponse of a large size quantum dot within time-dependent current density 
functional theory. Results for magnetic fields corresponding to integer filling 
factors are reported, as well as a comparison with the longitudinal dipole spin 
■ response. As in the two dimensional electron gas, the spin response at high 

spin magnetization is dominated by a low energy transverse mode. 
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I. INTRODUCTION 



Resonant inelastic light scattering has become a very useful tool to study quantum dot 
(QD) excitationsEHi. It had been extensively employed before to study electron single- 
layers and quantum wires in semiconductor heterostructures (see for example Refs. HH, 
and currently it is also being applied to study electron multi-layersBfi We refer the reader 
to Refs. for a review of this experimental technique. 

A major advantage of resonant inelastic light scattering over far infrared (FIR) optical ab- 
sorption techniques used in the pastBB is that the former allows to disentangle and identify, 
using polarization selection rules in the backwards geometry^, charge density (CDE), spin 
density (SDE) and single-particle (SPE) electron excitations in the same sample, whereas 
FIR absorption is only sensitive to charge density excitations. In inelastic light scattering, 
when the polarizations of the incoming and scattered photon are parallel (polarized geom- 
etry) CDE's are observed, whereas when the polarizations are perpendicular (depolarized 
geometry) SDE's dominate the spectrum. This is due to the structure of the scattering cross 
section, which besides the charge or spin electronic strength function contains the scalar or 
vector product of the photon polarizations, respectively^a. Yet, CDE's are seen with some 
intensity in the depolarized spectrumB. SPE's are mostly detected under conditions of ex- 
treme resonance, whereas CDE's and SDE's can be observed at incident photon energies far 
above the effective band gap!. This helps disentangle SDE's from SPE's, which are in the 
same energy range at small spin magnetizations. 

For a QD in the xy plane submitted to a static magnetic field B in the z direction, 
SDE's may involve electronic spin-flips or not. The later excitations are referred to as 
longitudinal SDE's, and the former as transverse SDE's. This means that SDE's are caused 
by one-electron excitation operators of the kind 

^ = E/(^K longitudinal (AS Z = 0) 

i 

(1) 

F± = ^2 fifi) a± transverse (AS Z = ±1), 

where S z is the z component of the spin of the dot, and a z is the z component of the Pauli 
matrix vector. <j± = (a x ±ta y )/2 are the spin-flip operators 

o+|T> = o--U> = 0; o+UHIT); o-_|T) = U). (2) 

Due to rotational invariance in spin space, longitudinal and transverse SDE's are degen- 
erated at B = if S z = in the ground state (gs). When the magnetic field is not zero, 
rotational invariance is broken by the Zeeman term. Then if the gs is almost paramagnetic, 
i.e. S z ~ as it occurs at even filling factors u, the SDE's are expected to split in a simple 
way: 

uX = uj lon ±g*fi B B , (3) 

where the superscript lon(tr) indicates the longitudinal (transverse) character of the mode, 
and the subscript ± corresponds to the two possible spin flip transitions, g* is the effective 
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gyromagnetic factor and \xb is the Bohr magneton. A thorough discussion of longitudinal 
dipole SDE's in quantum dots within time- dependent local-spin density theory (TDLSDT) 
has been presented in Refs. [3~5|JT6| . If the gs has a large S z , i.e. a large spin magnetization 
(ferromagnetic gs), as it happens at odd filling factors is, Eq. (^) no longer holds and longi- 
tudinal and transverse SDE's display dramatic differences arising from the spin dependence 
of the electron-hole (e-h) vertex corrections. 

To study transverse dipole spin modes we have resorted to current density functional 
theory (CDFT)ffl together with its time-dependent version (TDCDFT)§. For the physics 
we aim to describe here, local-spin density theory (LSDT) and CDFT sensibly yield the same 
results. CDFT is expected to be more reliable than LSDT at high magnetic fields§, but the 
residual interaction in the longitudinal spin channel is very cumbersome to work out within 
TDCDFT. For this reason, whenever we have had to work out the longitudinal spin response 
for the sake of comparison with the transverse one, we have resorted to TDLSDT. We will 
see that the spin dipole response is dominated by the transverse component, especially for 
filling factors smaller than 2. In the following we shall use the terms density functional 
theory (DFT) and time-dependent density functional theory (TDDFT) when the statement 
applies to either version of the general method. 

The spin and density response of the two dimensional electron gas (2dEG) has been 
thoroughly studied by Kallin and Halperin0 and MacDonaldH. As these authors, we have 
mainly addressed the response of quantum dots corresponding to integer filling factor gs's. 
We have considered only three cases whose ground state is the finite size analogue of a 
partially filled 2dEG configuration between v = 1 and 2. One should regard these results 
as qualitative, since the corresponding ground states are believed to be very complicated, 
strongly correlated ones, and the use of TDDFT to describe them may be questioned. 
Nonetheless, these results display the gross trends of the excitation spectrum, and for that 
reason we have considered them here. 



II. RESULTS 

We have taken of study a circularly symmetric QD of radius Rdot ~ 164 nm 

made of N = 210 electrons in a GaAs-AlGaAs heterostructure. The gs of this dot has been 



described in detail in Ref. £L8[ Throughout this work, we shall be using effective atomic 
units, whose definition and value for GaAs can be found for example in that reference. 

A major advantage of considering a rather large QD for the present discussion is that 
several integer filling factor gs configurations can be identified as a function of B. This 
allows to discuss the influence of the paramagnetic or ferromagnetic character of the gs on 
the excitation spectrum. The large number of relatively close single particle (sp) levels that 
one has to handle to obtain the gs and strength functions poses some technical problems, 
and is the token to pay for microscopically study large QD. In this work we have used a 
small temperature T < 0.1 K. This facilitates the calculations while it does not introduce 
any appreciable thermal effect in the results. 

In the following we limit the analysis to the most interesting dipole mode, for which the 
f(fi) in Eq. (|I|) is x,. Furthermore, to take advantage of the imposed circular symmetry, 
we have considered as dipole operators the following ones: 
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D±i fi = \ y £r j ^ai (AS z = 0) 
j 



2 i 



(4) 



as well as the combinations: 

5>V J Z = D +w + D_ li0 (AS 2 = 0) 

3 

(5) 

53^oi = D +1 , ± + D_ 1 , ± (AS Z = ±1) , 
j 

where the subscript ±1 represents the orbital angular momentum carried by the excita- 
tion, each one corresponding to a different kind of left- or right- circularly polarized light. 
Within TDDFT, the external probe represented by these operators has also a harmonic time 
dependence of frequency to. 



The method we have used to obtain the spin response has been described in Ref. IS for 



the transverse case, and in Refs. |15| , |16| for the longitudinal one, so we do not need to repeat it 
here. We want to recall that the mean field entering the Kohn-Sham (KS) equations changes 
due to the dynamical spin magnetization induced by the external field. Within TDDFT, it 
is assumed that electrons respond as free particles to the sum of the external plus induced 
field, which leads to a Dyson-type integral equation for the correlation function x 

X(r, r', lo) = x°(r, r u) + / dr x dr 2 x°(r, fi, u)V eh {r l: f 2 )x(r 2 , r u) , (6) 

where x° is the free electron correlation function and V e h is the residual e-h interaction. 
From x it is possible to obtain the induced spin magnetization corresponding to any of the 
above excitation operators and then determine the dynamical polarizability function a(u). 
Finally, the strength function S(u) is obtained from the imaginary part of the dynamical 
polarizability function as S(u) = -lm[a(uj)]. This procedure has been explicited in Ref. |19| 
for the -D+i,- operator. All these functions can be labeled according to the AL Z = ±1 and 
A5* z = 0, ±1 changes induced in the excitation process. In the following, we shall analyze 
the strength functions corresponding to the operators in Eq. (|5|), which we call xcr + , 
xa z and xa x = x(a + + cr_). 

Some characteristics of the strength functions are easily understood in terms of the 
uncorrelated e-h excitations (SPE's) used to build the corresponding x ^)? whose basic 
ingredients are the KS single particle energies e n i a and wave functions (p n ia(r) obtained 
from the solution of the KS equations which within LSDT read@ 

-iv 2 + 1 -u c £ z + 1 -u J y + V + (r) 



+ V H + V xc + (W xc + ^g*fi B B)a 2 



(fnia(r) = e n i a ip n i a (r) , (7) 



where V + (r) is the confining potential and V H = J dfp(f')/\f — f'\ is the Hartree po- 
tential. V xc = d£ xc (p,m)/dp\ gs and W xc = d£ xc (p,m)/dm\ gs are respectively, the spin- 
independent and spin-dependent exchange-correlation sp potentials obtained from the 
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exchange-correlation energy density S xc (p,m), and p and m are the electron and spin mag- 
netization densities p(f) = p^(r) + p^(f), m(f) = p\r) — p^(f). £ xc (p,m) has been 
constructed as indicated in Ref. [16|. Within CDFT, the KS equations are more cumbersome 
to write down0. Still, they have a similar structure, and what one has to keep in mind for 
the discussions that follow is the existence of an exchange-correlation potential which plays 
the role of V xc , and another exchange-correlation potential which plays the role of W xc in 
the above equation. 

One should be aware that the residual e-h interaction may change drastically the picture 
extracted from the uncorrelated e-h excitations. Yet, as a useful guide for the discussion, 
we collect in Fig. [L] the sp energies as a function of the angular momentum I for v — 8 to 
1, corresponding to B values from 1.29 to 10.28 T. In the bulk of the dot (small / values), 
these sp energies bear the band characteristics of the 2dEG, having similar e n i„ the sp states 
that have the same s z and yield the same Landau level index M = n + (|/| — l)/2, where 
n is the sp radial quantum number. For this dot, the filling factor v does represent the 
number of occupied Landau bands, each one labeled as (M, j) or (M, j). In Fig. [I], upright 
full triangles represent s z = j sp states, and downright, empty triangles, s z = j sp states. 
The horizontal lines represent the electron chemical potential. All the occupied bands are 
shown, but to build the valence space of sp states for calculating the correlation functions 
we have usually considered more empty states than those shown in the figure. 

It is worth to notice the small energy difference between the (M, f ) and (M, J.) bands 
for even v values AE^ ~ \g*psB\, which has its origin in the Zeeman term, as compared 
to the large energy difference between the same bands if v is odd, even if the applied B is 
relatively small; compare for instance the v = 7 and 6 cases. That difference in ferromagnetic 
gs's mostly comes from the spin-dependent exchange-correlation potential AE^ ~ 2|W^ 3:c |, 
which is zero or very small in paramagnetic gs's, and sizeable in ferromagnetic gs's, largely 
overcoming the Zeeman energy. Hartree-Fock (HF)H and LSDT or CDFTE3 yield such large 
gaps, whereas the Hartree approximation does not. The role of electron-electron interactions 



in producing these gaps when v is odd has been stressed in Ref. 24. We are going to see 



that the effect of the gap is paramount on the transverse spin response. 



A. Strength functions 

Figure || shows the strength function S{uj) (solid lines) corresponding to xa x . For the 
sake of clarity, we have decomposed xa x into its xa + and xo~_ components, which are drawn 
in Figs and £| respectively. The associated free responses, i.e. SPE's, are represented by 
dashed lines. The functions are given in effective atomic units, and the frequencies in meV. 

Let us first comment on the results corresponding to paramagnetic gs's in which both 
spin up and down sp states have a tendency to be equally populated yielding a small S z 
value. As a consequence, the attractive e-h residual interaction is weak, and SDE's and 
SPE's lie in the same energy range. In both xa + and x<7_ components, the strength displays 
a high energy structure with a frequency always close to the free strength (see also Fig. ^|), 
and a low energy structure. For paramagnetic gs's, the low energy excitation is a transverse 
spin edge mode built from e-h pairs near the Fermi level. These pairs can be easily identified 
in Fig. 1(a) and (b), as they are at the intersection of the chemical potential with the Landau 
bands. The sp band structure also explains why the edge mode is more fragmented at low 
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magnetic fields. For example, at v = 8 four e-h pairs, each one involving quite different sp 
orbital angular momenta from the other pairs, are contributing to the xcr_ strength, whereas 
only one pair is contributing at v = 2. These pairs are weakly correlated among them and 
the result is an edge mode fourfold fragmented at v = 8, threefold fragmented at v = 6, 
and so on. This nicely corresponds to the number of crossings of the Fermi level with the 
(Af, U) bands in Fig. |. 

In the case of xa + , the edge mode is less fragmented because some spin-flip e-h transitions 
having AS Z = +1 are Pauli blocked by our arbitrary election of B in the positive direction 
of the z axis (we recall that one has to have AL Z = ±1, which cannot always be fulfilled 
simultaneously with the spin and edge conditions). The lacking of the edge state in the v 
= 4 case is due to the particular sp structure around the Fermi level at B = 2.57 T. This 
accidental fact has no relevance for the general discussion. 

As anticipated, the e-h residual interaction produces a dramatic effect when the gs is fer- 
romagnetic. Even for moderately intense magnetic fields, the transverse spin mode emerges 
in these gs's as a very collective, undamped excitation whose energy we will see in the next 
Subsection depends little on the size of the dot provided it is of the present size or larger. 
This is somehow the analogous of Kohn's theorem^! for charge density excitations in a QD 
parabolically confined, and it has the same physical origin, namely, the exact (or nearly 
exact in the spin case) cancelation between the bare and induced e-h interactions. 

It can be seen from Fig. ^| that when the dot is fully polarized at v — 1, xo~ + no 
longer excites it because of Pauli blocking. Even at v = 3 {B = 3.43 T), its strength 
is very small. On the contrary, the excitation produced by icr_ displayed in Fig. f| is 
appreciably redshifted from the free response. The difference between both situations reveals 
the strength of vertex corrections arising from exchange-correlation terms of the electron- 
electron interaction, which within TDDFT00 are the only ones contributing to dress the 
free e-h vertex in the spin channel. This effect is more sizeable at v = 1 when the system is 
fully polarized. 

The low energy peak excited by xa_ is taking almost all the dipole strength, as we will 
show later. Low energy SDE's in ferromagnetic gs's caused either by xa + or xa_ are not 
edge, but bulk spin excitations. Again, Fig. [I] helps understanding this. For these odd v 
gs's, the Fermi level is between the (M max , j) and (M max , j) Landau bands, the former being 
occupied and the later empty. Although finite size effects distort the bands at the edge of 
the dot which is formed by sp states with high I values, it is clear that low energy spin- 
flip transitions involve sp states whose energy difference is precisely the energy difference 
between the (M max , f j) bands. These excitations also occur in the bulk (2dEG). The role 
of the residual interaction is clearly visible in Fig. [| comparing the free and TDCDFT 
strength functions at u—5 and at v — 3, for example. One sees that the SPE energies have 
nothing to do with the low SDE energy Thus, any free e-h model will be of little help to 
quantitatively analyze SDE's in partially polarized QD, the situation worsening the higher 
the polarization. 

We present in Figure [5] a more detailed picture of the spin excitation in the transverse 
channel at high B, showing the strength function of xa x in the 2 > v > 1 region. As in 
previous figures, the dashed lines represent the free response. In this range of filling factors, 
which corresponds to 5.14 T < B < 10.28 T, we have foundEl that the 25*2 value steadily 
increases^ from zero to 210, so the xa x strength is essentially that of xcr_ already discussed. 
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The interesting new feature in Fig. |5| is the structure of the high frequency peaks. They are 
two orders of magnitude less intense than the low energy ones, which thus exhaust most of 
the strength. Of the high energy peaks, the higher ones are excited by xcr_, and the lower 
ones by xa + (obviously, high energy transitions caused by xa + are blocked only when the 
system is fully polarized). It can be seen that these high energy peaks are little collective, 
as SPE's and SDE's are quite similar, and also that the centroid of the xa + and xo- peaks 
roughly follows the same evolution with B as the cyclotron frequency u c = eB/mc does. 
The value of u c is indicated in Fig. |5| by vertical arrows. 

When both high energy peaks are clearly visible in the strength, as for example at B 
= 7 T, their splitting is a quantitative measure of the spin-dependent exchange-correlation 
gap W xc , and its measurement may be the spectroscopic complement to experimental gap 
determinations based on the temperature dependence of the conductivity!?!!. This is so 
because W xc + g* /2 is directly related to the energy difference of the (M, f) and (M, {) 
bands around the Fermi level, see Eq. (|^). As the sets of "parallel" bands (M + 1, ][) and 
(M, f|) have the same "width" as Fig. p] indicates, the splitting between the high energy 
peaks would be twice the energy difference between the (M max , f) and (M max , J.) bands. 
This is nicely confirmed in the B = 7 and 9 T cases, for which an explicit comparison is 
possible (see Fig. 5 of Ref. [T8| ). We recall that this comparison is meaningful because of the 
weak effect of the residual interaction on the high energy peaks. 

The longitudinal high energy peak for these configurations lies@ at u lon ~ uj c (see also 
Fig. |7|), and an expression similar to Eq. (|3]) can be written: 

u tr _ Jon ±2 W XC . (8) 

The validity of Eqs. (0) and (|8|) is a consequence of the weakness of the residual e-h 
interaction when they hold (notice that g* and W xc are negative). 

Figure [3] shows the longitudinal (xa z ) and transverse {xa x ) dipole spin strength functions 
as dotted and solid lines, respectively. It can be seen that in ferromagnetic states, the 
strength is dominated by the transverse modes. One should also notice that for v > 2, i.e. 
low B, apart from some fine structure the longitudinal and transverse spin responses have 
their main peaks at quite similar energies. 

The paramagnetic configuration at v — 2 (B = 5.14 T) shows the interesting situation in 
which a zero spin gs sustains the simple result anticipated at the introduction: the transverse 
modes are just shifted by the Zeeman energy from the longitudinal ones. Interestingly too, 
the low energy xa + mode has almost collapsed. This suggests the presence of instabilities 
in this particular transverse channel, similar to those found in the longitudinal spin channel 
in the 2 > v > 1 regional. 

When the system is fully polarized, the longitudinal spin and charge density strengths 
trivially coincide and of the spin modes, only the transverse one has a senselll'@lll. The 
connection between the peak and Zeeman energies will be discussed in the next Subsection. 

The energies of the more intense, high energy peaks appearing in the xa + and the x<j- 
strength functions are shown in Fig. ^ as a function of B. The cyclotron frequency is also 
represented. Solid symbols correspond to even filling factor values, and open symbols to odd 
filling factor values from v — 10 (B — 1.03 T) to v — 1 {B — 10.28 T). Also drawn are the 
values of the high energy longitudinal spin peakslll (crosses). The triangles and diamonds 
correspond to the B = 7, 8 and 9 T peaks in the right panels of Fig. |5|. 
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In agreement with the preceding discussion, it is seen that for even filling factors the 
energies of the transverse SDE are ~ uo c ± g*fisB, thus close to u c , whereas for odd filling 
factors they are well apart from u c by the large spin-dependent potential W xc . In all cases, 
these peaks correspond to bulk modes involving interband e-h transitions made of sp states 
each one belonging to a Landau band with different index M. 

Similarly, Fig. |8| collects the energies of the more intense, low energy peaks appearing 
in the xa + (top panel) and of xo^ (bottom panel) strength functions. Solid and open 
symbols have the same meaning as in Fig. [7|. The Zeeman energy E z = —g*iisB > 
is also represented. To emphasize the energy staggering, consecutive v points have been 
connected by a thin line. As we have already discussed, these modes are spin edge modes 
for paramagnetic gs's, and spin bulk modes for ferromagnetic gs's. 

The SDE corresponding to xcr_ is the more interesting one. It is the only transverse spin 
mode that appears at high magnetizations, since the one generated by xa + is Pauli blocked. 
For ferromagnetic gs's, this is an undamped excitation since it is well apart from the SPE's 
(see Fig. |j). Notice that the transverse SDE energy approaches the Zeeman energy as B 
increases. At full polarization [y = 1), the energy is close to E z , but not equal to it. We 
shall come back to this point in the following. 



B. Sum rules 

Further insight onto the strength functions can be obtained from the evaluation of sum 
rules, which are their energy moments. Some of these moments are easy to obtain, model 
independent quantities. We are interested here in sum rules for non-hermitian excitation 
operators as the F± ones defined in Eq. ([!]). These sum rules have been extensively discussed 
in Ref. 



We consider the usual Pauli Hamiltonian H describing an N electron QD submitted to a 
constant magnetic field in the z direction (see for example Ref. JZ^) and define the following 
sum rules (SR): 

So = S -S+ = J2 |(n|F_|0)| 2 -E \(n\F + \0}\ 2 

n n 

(9) 

S 1 = Sr + St = $> n0 |(n|F_|0)| 2 +$> n0 |(n|F + |0}| 2 , 

n n 

where |0) is the gs of the system and |n) is an excited state with excitation energy u n Q. 
Using closure, it is easy to check that 

5 = (0|[F +) F_]|0> 

(10) 

5 1 = (0\[F + ,[H,F_]\0) 
An explicit evaluation of So and Si yields 

5 = J df\f(f)\ 2 m (f) 

(11) 

51 = l -J df\Vf(r)\ 2 Po(r)-9^BB J dr\f(r)\ 2 m (r) , 



S 



where po(r) is the gs density of the dot, and m (r ) = p (f) — Po(r) is the gs local spin 
magnetization. These sum rules are relating properties of the exact spectrum of the Pauli 
Hamiltonian to properties of its exact gs. 

So and Si are also fulfilled within TDDFT. Indeed, it can be proved using the techniques 
discussed in Ref. that the TDDFT spectrum is such that one can obtain them with 
TDDFT accuracy from Eqs. fllTD using the KS gs. Not all approximation schemes fulfill 
these sum rules. Independent particle spectra such as those obtained from HF, KS or Hartree 
approximations violate Si. This means that one would not obtain the second of Eqs. flllD 
using in the second of Eqs. (|It]) the corresponding one-body Hamiltonian. For example, the 
KS spectrum we use to build the SPE's and x° correlation function yields 

S^ s = ^ J df\Wf(n\ 2 Po(n ~ J dr[g*p B B + 2W xc (f)}\f(f)\ 2 m (f) . (12) 

Within TDDFT, the effect of the e-h induced interaction is crucial in restoring Si: it exactly 



cancels the 2W XC contribution in Eq. (12). We want to point out in passing that in the 
density channel, the induced interaction is also responsible for the fulfillment of Kohn's 
theorem which would be violated otherwise. 

When the dot has a large gs spin magnetization, the xa + term in Eq. (|9|) contributes 
very little to So and Si. Thus, one is left with the xa_ contribution, which is concentrated 
in a narrow energy region (see right panels in Fig. [|). Under these conditions, it is a fair 
approximation to identify to = Si/S with the mean energy of the peak displayed in the icr_ 
strength function. To proceed further, let us consider that the dot is fully polarized, i.e., 
^o(^) = Po{r)- Taking f{f) = x, one gets 

50 = S- = !(r 2 > 

(13) 

C <?- N N * W 2\ 

51 = S l = — - —g p B B{r ) , 
where (r 2 ) is the root mean square radius of the dot 

(r 2 ) = — Jdfr 2 p {f) . (14) 

Thus 

o = ^ = -g*l* B B + -±r. (15) 
S (r*) 

At v = 1, taking for poif) that of the maximum density droplet (MDD) stateS, which is a 
good approximation for large dots, we have 

s _| = _ 9 . MbB + _*_. (M) 

Expressions (p~5|) or (p^6|) are the SR estimates of the transverse dipole SDE at large spin 
magnetization. Using again the MDD (r 2 ) value, we get Sq ~ N(N + l)/2u c . This shows 
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that for high spin magnetizations, the squared xcr_ matrix elements in the SDE's are a 
factor ~ N stronger than the squared x matrix elements in the CDE's, whose strength is 

For GaAs, u c /\g*fi B B\ = 2/\g*m*\ ~ 68. We thus have Q ~ 0.35 meV for the JV = 210 
dot at v — 1, whereas the peak energy is E pea k ~ 0.27 meV and E z ~ 0.26 meV. In the 
limit of a very large dot, Q becomes E z . This is the analogous of the 2dEG result§@ for 
the spin- wave dispersion relation at q = (Larmor's theorem). 

So far, we have discussed the response to the dipole L = 1 operator. It is straightforward 
to consider the response to a general L mode operator of the kind f(r) ~ r L e ±tLe . These 
fields are relevant to study spin and charge density modes with well defined angular mo- 
mentum. Moreover, in Raman experiments one may face a situation in which the in-plane 
transferred linear momentum q is small enough, so that qRdot << 1 and the plane wave 
operator e iqr involved in the excitation process can be expanded into r L e ±lL0 multipoles, 
each of them probing a well defined angular momentum CDE or SDE0. When q cannot be 
considered as small!, the TDDFT response may still be worked out, fixing q and adding the 
responses to f(r) = JL{qr)e ±tLe , because e tqr = Y1,l iL Jiiqr) e lLd . The number of terms 
in the expansion may be large depending on the q value, but in principle the method is of 
direct applicability. 

We finish this Subsection discussing the SR for f(r) = e iqr in the more interesting case 
of full polarization. Since \f(f )| 2 = 1 and |V/(r )| 2 = q 2 , we get from Eqs. ( p]) 

5 = N 

(17) 

51 = ^q 2 N - g*fi B BN . 

The same equations hold for the 2dEG substituting N by the electron density. In either 

case, 

u = -g^ B B+ l -q 2 (18) 



We want to emphasize that Eqs. fllT]) and (|18D are exact, model independent and valid 
for any q value. Only the applicability of the Pauli Hamiltonian to describe this physical 
situation has been taken for granted. 

At first glance, Eq. (|TSj) is in contradiction with the spin-wave dispersion relation of 
Ref. pHP2"|, whose q dependent term has an electron-electron energy dependence instead of 
the 1/2 factor. The difference stems from the sp valence space, which is different in both 
calculations. Indeed, the sum rule result Eq. (|TjI) takes into account all possible intraband 
and interband sp excitations induced by the operator e iqr ' cr_. Thus, uo is an average of 
the low and high energy peaks. In Refs. |21|J22] , the valence space was restricted to the filled 
(0, t) and empty (0, {) bands to specifically address the low energy mode. Their result can be 
exactly recovered in the SR approach if one uses the same valence space and accordingly, the 
projection of e iqr ct_ onto the M = spaced. This is the so-called single mode approximation 
(SMA)il, equivalent to the approach of Ref. P]],P2| in the fully polarized case. 



It is easy to seize the effect of the M > bands on the low energy collective mode within 
TDCDFT. It suffices to compute the response limiting the sp valence space to the M = 
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bands. In this case, only the low energy peak appears in the xcr_ strength function, and its 
energy ~ 0.30 meV is denoted by a cross in the bottom panel of Fig. |8|. The inclusion of 
the high energy bands changes the energy of the transverse SDE in ~ 10 %. We have also 
determined that the high energy peaks in the xcr_ strength contribute around 20 % to Si. 



III. SUMMARY 

In this work we have thoroughly studied the transverse dipole spin response in quantum 
dots. Together with our previous worksEl'Ei@, they provide a detailed account of the appli- 
cability of time-dependent density functional theory to the description of the dipole response 
of QD in the charge and spin channels. Although microscopic descriptions of similar 
plexity exist since some time ago for the charge density modes in medium size QDt^ 
the spin density modes in QD had not been previously addressed. Besides, we have been 
able to applythe theory to rather large dots, similar in size to those investigated in present 
experimentalHlllra 

Among SDE's, the transverse ones are especially relevant; in the longitudinal channel, the 
residual interaction is fairly weak, and the SDE's are Landau damped as they are close to the 
SPE's (actually, the same happens with the transverse modes at low B). In the transverse 
channel, when the dot has a sizeable magnetization the position of SDE's is shifted away 
from the SPE's by exchange-correlation vertex corrections arising from electron-electron 
interactions. As a consequence, a very collective, dispersionless SDE emerges. At large spin 
magnetization, Pauli's principle plays a prominent role, blocking the xa + component of the 
transverse strength function, which becomes simpler. 

The possibility of carrying the calculations in a large dot displaying several integer quan- 
tum Hall gs's has permitted us to disclose the sensitivity of the transverse response to the 
applied B, which appears as a strong energy oscillation with B, and a nearly collapsed low 
energy xa + mode at the v = 2 paramagnetic gs. The energy oscillations are also consequence 
of the different strength of vertex corrections in ferromagnetic and paramagnetic gs's. 

We have also derived two model independent sum rules that, on the one hand, can be 
used as a control for the analysis of the experimental data, and on the other hand, have 
allowed us to relate our results at full magnetization to previous works on 2dEG. 
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FIGURES 



FIG. 1. (a): Single-particle energies as a function of orbital angular momentum I for v = 8 to 
5. The horizontal solid lines represent the electron chemical potential. The full, upright triangles 
represent (M, f) bands, and the empty, downright triangles represent (M, 1) bands, (b): Same as 
(a) for z/ = 4 to 1. 

FIG. 2. Strength function corresponding to xa x (solid lines). The dashed lines represent the 
free strength function. S(u>) is in effective atomic units divided by 10 5 . 

FIG. 3. Strength function corresponding to xcr+ (solid lines). The dashed lines represent the 
free strength function. S(u>) is in effective atomic units divided by 10 5 . 

FIG. 4. Strength function corresponding to xo- (solid lines). The dashed lines represent the 
free strength function. S(u>) is in effective atomic units divided by 10 5 . 

FIG. 5. Same as Fig. 2 for 2 > u > 1. The vertical arrows indicate the value of uo c . 

FIG. 6. Strength function corresponding to xo x (solid line) and xa z (dotted line). The strengths 
are in effective atomic units divided by 10 5 in the transverse case, and by 5T0 4 in the longitudinal 
case to make it easier to distinguish them. 

FIG. 7. Energies of the more intense high energy peaks excited by xa + (squares) and xo- 
(circles) as a function of B. Solid symbols correspond to even filling factor values, and empty 
symbols to odd filling factor values from v = 10 (B = 1.03 T) to v = 1 (B = 10.28 T). Also drawn 
are the energies of the high energy longitudinal spin peaks (crosses). The triangles and diamonds 
correspond to B = 7, 8 and 9 T. 

FIG. 8. Bottom panel: Energies of the more intense low energy peaks excited by xd_ as a 
function of B for the same configurations as in Fig. 7. Also shown is the Zeeman energy (dashed 
line). The cross at u = 1 is the SDE value obtained from TDCDFT when the sp valence space is 
limited to the M = bands. Top panel: Same as bottom panel for xa + . Some configurations are 
absent due to Pauli blocking. 
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